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Abstract. The Baxter-Bazhanov-Stroganov model (also known as the r'^^ model) 
has attracted much interest because it provides a tool for solving the integrablc chiral 
Zjv-Potts model. It can be formulated as a face spin model or via cyclic L-operators. 
Using the latter formulation and the Sklyanin-Kharchev-Lebedev approach, we give 
the explicit derivation of the eigenvectors of the component Bn{X) of the monodromy 
matrix for the fully inhomogeneous chain of finite length. For the periodic chain we 
obtain the Baxter T-Q-equations via separation of variables. The functional relations 
for the transfer matrices of the r^^^^ model guarantee non-trivial solutions to the Baxter 
equations. For the N = 2 case, which is free fermion point of a generalized Ising model, 
the Baxter equations are solved explicitly. 
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1. Introduction 

The aim of this paper is the exphcit construction of eigenvectors of the transfer- 
matrix for the finite-size inhomogenous periodic Baxter-Bazhanov-Stroganov model 
(BBS model) also known as the r(2)-model P 121 El ij. This is a iV-state spin lattice 
model, intimately related to the integrable chiral Potts model. The connection between 
the 6-vertex model, the BBS model and the chiral Potts model gives the possibility to 
formulate a system of functional relations jSlIll for the transfer matrices of these models. 
Solving these systems is the basic method for calculating the eigenvalues of the transfer 
matrix of the chiral Potts model [S] , and under some analyticity assumptions, to derive 
the free energy of this model [S] and its order parameter [Hj. 

In general for the BBS model there is no Bethe pseudovacuum state and so the 
algebraic Bethe ansatz cannot be used. Therefore, in order to achieve our goal, we shall 
use the formulation of the BBS model in terms of cyclic L-operators first introduced by 
Korepanov [Zj and Bazhanov-Stroganov P] and adapt the Sklyanin-Kharchev-Lebedev 
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method of separation of variables (SoV) [HI El UHl HI] for solving the BBS eigenvector 
problem. The fusion equations will provide the existence of solutions to the Baxter 
equations. 

The paper is organized as follows. After defining the BBS model as a statistical 
face spin model, we give the vertex formulation of the model in terms of a cyclic L- 
operator and conclude the Introduction explaining the two basic steps involved in the 
SoV method. Section 2 deals with the solution of the auxiliary problem, leaving for 
Section 3 the lengthy inductive proof of the main formula. The short Section 4 derives 
the action of the diagonal component D of the L-operator on the auxiliary eigenstates. 
Then in Section 5 we come to the periodic model and in deriving the Baxter T-Q- 
equations we show the role of the fusion equations for solving these Baxter equations. 
In Section 6 we apply these results to the homogenous N = 2 case and calculate the 
eigenvalues and eigenvectors of the homogenous N = 2 BBS model which is the free 
fermion point of a generalized Ising model. Section 7 gives our conclusions. In an 
Appendix we show a strong simplification occurring if the BBS model is homogenous. 

1.1. The BBS model 

Following the notation of a recent paper of Baxter we define the BBS model as a 
statistical model of short-range interacting spins placed at the vertices of a rectangular 
lattice. We label the spin variables s^:,^ by a pair {x,y) of integers: x = 1, . . . ,n, and 
y = 1, . . . , m. Each spin variable Sx,y takes N values (A^ > 2): 0, 1, . . ., A^ — 1. The 
model shall have Zjy-symmetry and we may extend the range of the spins Sx^y to all 
integers identifying two values if their difference is a multiple of A^. The model has a 
chiral restriction on the values of vertically neighboring spins: 

Sx,y - Sx,y+i = or 1 mod A^ . (1) 

In the following we will consider the spin variables on two adjacent rows: {k, I) and 
(fc, / + !), where I is fixed and k = 1, . . . ,n. Let us denote Skj = 7^ and s^j+i = 7^. The 
model depends on the parameters tq and a'^, ,b'i. ,c'j^ ,d'i., a'l ,b'l ,c'l ,d'l, k = 1,2, ... ,n. 
Each square plaquette of the lattice has the Boltzmann weight (see FiglH) 

1 

l^.(7fc-i,7.;7Li,7D= E uj^''-'^'''^-''-'\-ujtqy>'-^'^-^--^x (2) 
x^fc-i(7fc-i - 7fc-D"^fc-i)^fc (7fc - 7fc,"^fc-i), 

where to = e'^'^^^^ , and 

i^^(0,0) = l, F^(0,l) = -^t,^, F^(1,0) = ^, Fl{l,l) = -u:f, 

"k "fc "k 

and expressions for F^'(7fc — 7^, rrik-i) are obtained from ^^(7^ — 7^, nik) by substitutions: 

n' h' r' d' ^ n" h" r" d" 
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Figure 1. The triangle with vertices markes by the spin variables 7^-1, 7^;_i, mk-i 
correspond to the function F^_]^(7fc-i — 7fe_i,?7ifc_i) in the triangle 7fc, 7^, irik-i 
to Fl!{jk - 4,mfc_i). 

We will consider the periodic boundary condition: 7^+1 = 71, ^'n+i = where n 
is the number of sites on the lattice along the horizontal axis. The transfer-matrix of 
the periodic BBS model is A^" x A^" matrix with matrix elements 

n+l 

tn(7,7') = llWr{^k-l,Jk;iLl,l'k), (3) 
k=2 

labelled by the sets of spin variables 7 = {71, 72, ... , 7„} and 7' = {7J, 72, ... , 7^} of 
two neighbour rows. 

Considering m^, /c = 1, . . . , ra, in (j2I) as auxiliary spin variables which take the two 
values and 1, we can rewrite transfer-matrix Q in a vertex formulation associating a 
statistical weight not to the plaquettes but to vertices each of them relating four spins: 
nik-i, nik, 7fc, 7fc (see Fig. 1). Then the weight associated with the A;th vertex is 

4(t,;mfc_i,mfc;7fc,7^) = ^™.-i7;-™.7.(_^i^)7.-7^— .-1 x 

X ^kilk - Ik, mk^i)Fl^{-fk - Ik, mk). (4) 
and the transfer-matrix (jSj) can be rewritten as 

n+l 

tn(7,7')= ^ Ylh{tq;mk-i,mk;-fk,'yk)- (5) 

mi,...,m„ k=2 

1.2. The L- operator formulation o/BBS model 

For our construction of the BBS model eigenvectors we will use a description of this 
model as a quantum chain model as introduced in [7j and [3]. To the each site k of the 
quantum chain we associate the cyclic L-operator acting in a two-dimensional auxiliary 
space 
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Lk{\) 



Ufc(cfc-4vfc), XakCk + ^khdk/ Xk 



k = l,2,...,n. (6) 



At each site k we define ultra-local Weyl elements and obeying the commutation 
rules and normalization 



u 



N 



1. (7) 



In (jn}, A is the spectral parameter and we have five parameters x^, Ofc, 6^, c^, rf^ per 
site. At each site k we define a iV-dimensional linear space (quantum space) Vk with the 



basis |7)fc, 7 G T^n and natural scalar product a:(7'|7)a: 
Ufc and Vfc act by the formulas: 



5^/,^. In Vk the Weyl elements 



(8) 



^k\l)k = ^^^\l)k, Vfc|7)fc = |7 + 

The correspondence between the lattice BBS model and its quantum chain analog 
is established through the relation 



ikitg;mk-i,mk;'jk,7k) = k{7\Lk{X) 



|7')* 



(9) 



and the following connection between the parameters of these models 



A = —ujt. 



Xk 
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We extend the action of the operators u^., to V*^"^ = Vi ®V2 ® ■ ■ ■ ®Vn defining 
this action to be trivial in all Vs with s ^ k. The monodromy matrix for the quantum 
chain with n sites is defined as 

a,(A) /^„(A) 

The transfer-matrix is obtained taking the trace in the auxiliary space 



T„(A) = Li(A)L2(A) ■■■L„(A) 



(10) 



t„(A) = trT„(A) = A„(A) + D„(A). 

This quantum chain is integrable because the L-operators 
twisted 6-vertex i?-matrix at root of unity 

f \ -ujv 

uo{\-u 

v{l-uj 





R{X,u) 



\ 





X{l-uj) 
A-z/ 




(11) 

are intertwined by the 
\ 



A 



ujiy 



J 



i?(A,z/)L« (A)L 



'(2), 



Lf{u)Ll'>{X)R{X,u), 



-(1) 



(12) 



(13) 



where L'j}\X) = Lk{X)m, (A) = I®Lfc(A). Relation (0 leads to [t„(A), t„(/i)] = 
and so t„(A) is the generating function for the commuting set of non-local and non- 
hermitian Hamiltonians Hq, . . . , H„ of the model: 

+ H„_iA"~"^ 



t„(A) = Ho +HiA + 



HnX'' 



(14) 
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The lowest and highest Hamiltonians can be easily written explicitly in terms of the 
global Zjv-charge rotation operator V„ 



n 



Ho = 1 +V„]^^^; H„ = JJ^ OfcCfc + V„]^Xfc; V„ = viV2---Vn. (15) 

k=l ''^^ k=l k=l 

It also follows from the intertwining relation (jlHj) that Bn{X) is the generating function 
for another commuting set of operators hi, . . . , h„: 

[5„(A),5„(^)] = 0, 5„,(A) = hiA + h2A2 + ... + h„A". (16) 

The great interest in the BBS chain model is due to its relation to the integrable 
chiral Potts model. In [HIIZI it was observed that besides the intertwining relations (fT!^ . 
the L-operators © satisfy a second intertwining relation in the Weyl quantum space: 

Sc,ia2;/3i/32 (p, p\ g, q') ^f/' P) Ljl'J' (A; g, q') 

/Sift J 

= J^I^r^f (A;g,g')^-f (A;P,P') S,,^,;,,,,(p,p',g,g'), (17) 

/3i/32J 

if the parameters are chosen as 

Xfc = , afc = Xq^^, bk = , Ck =uJXg'^, dk = , (l»j 

f^Qkf^g'k f^qkf^g'k f^gkf^q'k 

where Xg^., i/g^, (analogously for Xqi^ etc.) satisfy the chiral Potts model constraints 
<+< = k(« + l); k< = l-kV--; k< = l-kV-; k^ + k'^ = 1. (19) 
Here k and k' are temperature parameters. In ()17|) we have written the spin matrix 
elements of the L-operators with the parametrization (|T8|l as L'^f {\,q,q') where 
j = 0, 1 are the components in the auxiliary space and greek indices q;,/? = 0,...,A^ — 1 
denote the components in the quantum space (jH)), suppressing the site index k. The 
matrix S turns out to be the product of four chiral Potts-Boltzmann weights jB] 

Sa^a2Ai32{p^p\ q, q') = Wpg>{ai - a2)Wp'g{p2 " Pi)Wpg{p2 - ai)Wp'g,{pi - aa). (20) 
In the parametrization (fTHj) of the BBS model there are various functional relations to 
the chiral Potts model transfer matrix which have been used to obtain explicit solutions 
for the chiral Potts eigenvalues |3]. Only further restricting the parameters to the 
"superintegrable chiral Potts" case: 

ak = Lo'^bk = Ck = dk = Xk = I (21) 

allows to solve the BBS model by algebraic Bethe ansatz, see e.g. ^SIEIE]- In this 
form Baxter in 1989 first obtained the BBS model as an "inverse" of the superintegrable 
chiral Potts model, see Equations (8. 13), (8. 14) of j2j. In this paper we shall follow P] 
in not using restrictions like (fT^ on the parameters >Ck, ak, bk, Ck, dk- We only shall 
exclude the superintegrable case (pT|l . 

It was shown in the paper [THj that the N = 2 BBS model is equivalent to the 
generalized Ising model at the free fermion point. The results of our paper permit to 
obtain the transfer-matrix eigenvectors for this model. Recently, the interesting paper 
PTj appeared, where these eigenvectors was constructed using the Grassmann functional 
integral. 
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1.3. Functional Bethe ansatz and SoV 

The construction of common eigenvectors of the set commuting integrals (|T4jl will be 
solved in two main steps, which generally can be formulated as follows: 

First, for the given quantum integrable chain type model one has to find an auxiliary 
integrable model such that: 1) the eigenvectors for original model can be expressed as 
a linear combination of the eigenvectors for the auxiliary model and 2) the coefficients 
of this decomposition should factorize into products of the single variable functions 
(phenomenon of "separation of variables"). 

Second, the auxiliary problem should be chosen in such a way that the construction 
of its eigenvectors is a simple iteration process: eigenvectors for the auxiliary model of 
size n have to be obtained from the eigenvectors for the model of size n — 1. 

An example of realization of the first step in case of the Toda chain model was 
proposed in the paper ^HJ- The auxiliary model for the periodic Toda chain was the open 
Toda chain. In the paper ^H] on the Toda example this approach has been formalized 
as "the functional Bethe ansatz" jH] . A complete realization of this step for the periodic 
Toda chain model can be found in the paper |S|. 

Regarding a recurrent procedure for eigenvectors of the auxiliary problem, probably 
the first reference to this possibility may be found in the series of lectures [20] • The 
main idea of this approach can be formulated as follows. Consider an integrable quantum 
chain model of size n. The monodromy matrix is the product of n L-operators. We 
decompose a system into two subsystem of the sizes ni and n2 such that n = ni + n2. 
Suppose we can solve the eigenvalue and eigenvector problems for the subsystems. In 
PU] E. Sklyanin claims that there is a relation between the eigenvectors of the original 
system and the eigenvectors of its smaller subsystems. For the open Toda chain model 
this Sklyanin approach has been realized in the paper [TUj . 

In our case of the BBS chain model, the auxiliary model is governed by the set of 
commuting integrals ()16|). So we first solve the problem of finding the eigenvectors for 
these integrals. Then we shall show that the eigenvectors for the operators ()14|) can 
be constructed as linear combinations of the eigenvectors of the set fllfij) . The multi- 
variable coefficients of this decomposition admit the separation of variables and can 
be written as products of single variable functions, each satisfying a Baxter difference 
equation. We shall obtain this Baxter equation for generic and solve it explicitly for 
N = 2 corresponding to the free fermion point of the generalized Ising model |ini. Note 
that the eigenvectors of the commuting set of operators which come from the generating 
polynomial An{X) {[An{X) , An{fi)] = 0) were found in paper pT| . 

2. Eigenvectors of Bn{X). 

2.1. Consequences of the RTT relations 

We start with the second step of the program described in Section 11.31 Following 
[HEnilin], we first construct eigenvectors of Bn{\). According to (|T6|) . any common 
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eigenvector of the commuting set of operators hi, . . . , h„ is eigenvector of Bn{\) and 
the eigenvalue is a polynomial in A. Factorizing this polynomial we get 

n-1 

5„(A) = AAo n - ^'^) ^ = (^0, Ai, . . . , A„_i), (22) 

k=l 

where we labelled the eigenvector \1/a by the normalizing factor Aq and the n — 1 non- 
vanishing zeros Ai, A2, . . • , A„_i of the eigenvalue polynomial. 

Now the intertwining relations (fT!?|l tell us that if is eigenvector of -B„(A), then 
by applying repeatedly the operators A{Xj) and -D(Afc) {j,k = 1, . . . ,n—l) and V„ ffT3j) . 
we can generate a whole set of A^" eigenvectors of i?„(A). 

The intertwining relations ()13|) give 

(A - iUfi)A^{X)BM = ^(A - /i)5„(/i)A„(A) + /i(l - iu)AMBniX) (23) 
(A - ujfi)DMBniX) = ^(A - fi)B^{X)DM + A(l - uj)D^{X)BM. (24) 
Fixing A = X^, k = 1, . . . ,n — 1 in ()23|) and acting by it on \E'a we obtain 

BM {An{Xk)^x) = f^Xoifi - io-'Xk)UsM^' - iM>^k)^x) ■ (25) 
This means that up to a constant 

A„(A,)M/A~^Ao,...,a.-lA„...,A„_i- (26) 

Similarly, from ()24|1 we can get 

^n(Afc)^A ~ ^^-iAo,...,..A„...,A„_i. (27) 

(The proportional factors will be obtained later in (fTOI) and fl90p ). Furthermore, acting 
by ()23|1 on \1'a and extracting coefficients of A"^^/i"' we have 

VA~^.-iAo,...,A„...,A„_i- (28) 

We see that the operators A„(A) and -D„(A) at the eigenvalue zeros A^ of -B„(A), 
together with the charge rotation operator V„ = V1V2 ■ ■ ■ v„ act as cyclic ladder 
operators on the eigenvectors of -B„(A). So the eigenvalues of -Bn(A) can be written 

n-1 

E„(A)^p„ = A r„,oo;-''-« J] + ^n,kOO~'"'') > (29) 

k=i 

where Vn^s, s = 0, . . . , n — 1 is a set of constants and we shall use the phases 

Pn = (Pn,0, • • • , Pn,n-l) ^ {Z^T (30) 

as new labels of the eigenvectors. The fact that the eigenvectors of the operator Bn{X) 
can be considered as dependent only on the integer phases of the roots 

An,fc = -rn,k^~^"''' (31) 

is a common property of the root of unity integrable models. The amplitudes r„ ^ of 
these roots are fixed by some "classical" procedure which will be described below. In 
some cases this procedure becomes an classical integrable system naturally incorporated 
into the quantum system (see [201 and references therein). 
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2.2. One- and two-site eigenvectors for the auxiliary problem 

We now start to solve the auxihary problem, which is to compute the eigenvectors 
of Bn{\) in the basis V*^""-*. We shall adapt the recursive procedure of Kharchev and 
Lebedev jTU] to the BBS chain model. 

In our root of unity case a very important role will be played by the cyclic function 
'W^p(7) |2H| which depends on a Z^v- variable 7 and on a point p = {x, y) restricted to 
the Fermat curve + y^ = 1. We define Wp{'~f) by the difference equation 

Wph) y N , N 



Wp(7 — 1) 1 — CJTX 

The Fermat curve restriction guarantees the cyclic property Wp{^ + N) = Wp{'y). The 
function ^^(7) is a root of unity analog of the g-gamma function. 

It is convenient to change the bases in the spaces Vk- Instead of |7)a;, 7 G Z^, we 
wiU use the vectors 

V'f = E - P^^^^ (33) 

which are eigenvectors of the upper off-diagonal matrix element Xu^^^a^. — fefoV^) of the 
operator L^. 

= XakY^ Wp^{-f - p)uj-^\-f)k - XbkY^ ti;p,(7-p)^"(^+^)|7 + l)fc 



X y^Wp^i^y-p) (ak- — ]uj '' + bk—uj' 



" \^), = XTuuj-'^i,f\ (34) 

\ III. J iiu 

In the first step we used (jH)) and to obtain the last line we used with y^ = b^/ak, 
Tfc = Xkttk- The Fermat curve restriction for p^ = {xk,yk) gives = — b^ . We see 
that if Tfc = (in particular, in the superintegrable case (j2II)) it leads to Xk = 0, yk = 1- 
In this case ()33p does not give a new basis in Vk- This is a reason why we exclude values 
of the parameters which lead to the degeneration of the cyclic function ^^(7). 

This sequence of operations applied in ()34|1 will be performed rather often in the 
following derivations. The application of shifts the spin index. This is compensated 
by the shift of the summation variable, which results in an opposite shift of the argument 
of Wp. This in turn is removed using ()32|) . 

The operator Vfc shifts the index of tp^''^: 

v^V'f = 5^ «^p,(7-p)l7 + l)fc= J2'^P^(^-p-^^\^)'^ = 4+i- (35) 

Using (jH^ for = 1 and comparing to (f^. we write one-site eigenvector as ^ '■= 
i'tifi- With Tifl = ri we have 

5i(A)^p,„ = A ri,o u;-P''° ^p, , , Ai(A)^p, = ^p, , + Axi^^, „+i . (36) 
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The construction of the two-site eigenvectors will show us the first step of 
the recursive method. In accordance with ()29|) we are looking for eigenvectors \E'p2, 
(p2 = (p2,o,P2,i) € ^Af X ^Af) of the two-site operator B2{\), which should satisfy 

52(A)^p, = A r^^o^-P^'^' (A + r2,i^-^^'^)^p, . (37) 

We suppose that ^ can be written as a hnear combinations of products of one-site 
eigenvectors 

^P. = E Q(Pi'P2|P2)^i>V^(?. (38) 

P1,P2&N 

Using and (jBll), the matrix Q can be calculated as follows: 
52(A)^P2 = (Ai(A) A u^^(a2 - &2V2) + 5i(A)(Aa2C2 + 62 c^2V2/x2)) ^ p^ = 



P2 



= E {<5(pi,P2|p2) (Ar2t^ + \'^a2C2riu + Q{pi - 1, p2|p2)A^xir2t^" 

Pl,P2 

+ g(pi,p2 - l|p2)^Ari^-^4^ii^® ^i?- (39) 

^2 J 

Comparing powers of the spectral parameter A in pHjl and in (jSZ|), together with pHjl . 
we get 

(r2,o^~^^'°-a2C2ria;-''i)Q(pi,P2|p2) = xir2a;-''^Q(pi - 1, P2IP2), (40) 

(r2,or2,i^-''^'°-''^-^ -r2aj-''2)Q(pi,p2|p2) = ^riu;-''^Q(pi, P2 - l|p2). (41) 
The difference equations pUj) and ()41|) have the solution 

^-(P2,0+P2,l-Pl){P2,0-P2) 

Q(P1,P2|P2) = -T ^ ■ -T, (42) 

^P2,o(P2,0 - Pi - l)Wp2(p2,0 + P2,l - P2 - 1) 

where p2,o = (a;2,o, 1/2,0), V2 = (^2,^2) and 

'^i 1^2 _ r2 _ &2 (^2 '^l 

X2,0 = 0.2 C2 , 1/2,0 = ><1 , 2:2 = , 1/2 = • (43) 

^^2,0 ^^2,0 ^^2, 0^^2,1 X2 r2,o?^2,l 

The parameters r2,o and r2,i are determined from the condition that the points p2,o and 
P2 belong to the Fermat curve. 

One can proceed this way to construct n-site eigenvectors of the auxiliary problem. 
In fact, in order to see the general structure emerging, one has to go to the four-site 
case. We shall not do this here, but rather in Section IHl we shall prove the general 
result by induction. This proof will use recursive relations between amplitudes r„fc, 
k = 0,1, . . . ,n — l formulated in the Subsection 12.41 Fermat parameters p = {x, y) of the 
cyclic functions Wp{p) appearing in our construction will depend on these amplitudes. 
Compatibility conditions between recursive relations for the amplitudes and Fermat 
curve equation x^ + = 1 can be formulated as a "classical" BBS chain model j29j . 
This model will be formulated in the next Subsection using an averaging procedure for 
the cyclic L-operators (jH)) given in 
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2. 3. Determination of the parameters r^ s 

Let us define the "classical" counterpart O(A^) of a quantum cyclic operator 0(A) 
using averaging procedure [Hj 

O(A^) = (0)(A^)= n.ez.OKA) (44) 

and apply this procedure to the entries of the quantum L-operator (jHl). Denote the 
result by Ck{X^) 

\{L,,) (Ln) ) \ cjr-< - eA'-of <f ; 

where e = (—1)^, and call it as the "classical" ^-operator of the classical BBS model. 
Accordingly, the classical monodromy for the m-site chain is 

r„(A-) = A(A-)A(A-).../:^.(A-) = I ^-^^2 f'"^.!!! 1 (46) 



Cm,(A^) T>m{)\^^ 

where the entries are polynomials of A^. By Proposition 1.5 from |14j, these polynomials 
coincide with averages (Am), (Bm), {Cm) and (Dm) of the entries of (|Tn|) . This 
proposition provides a tool for finding the iV-th powers of the amplitudes r^.^: applying 
(HH) to JSni) we obtain 

m—l 

Bm{X^) = (-6)-A^<o n (^"^ - ^ <s) ■ (47) 

s=l 

This relation together with and ()4(i|l allows to find ^ in terms of the parameters 
, , , and x^, k = 1, . . . , m. The problem of finding the amplitudes rm,s is 
reduced to the problem of solving a (m — l)-th degree algebraic relation. As shown in 
the Appendix, in the case of the homogeneous BBS chain model the problem is reduced 
to solving a quadratic equation only. The described procedure gives the amplitudes rm,s 
up to some roots of unity. In fact we can fix these phases arbitrarily because this leads 
just to relabeling of the eigenvectors. In what follows we suppose that we fixed some 
solution {rm,s} in terms of the parameters , b^ , , d^ and x^. 

Let us give a recursive description for Bm{X'^)- From ()46|). we immediately read off 
the recursion relations 

^™(A^) = (1 - ex^A^)^„_i(A^) + (c^ - rf^)S„_i(A^) , 

i3^(A^) = -eA^(a^ - 6^)An-i(A^) + (bX/xZ - eA^a^c^)i3™_i(A^) . (48) 



Combining these two relations we get 



^^(A^) = ' ^^(A^) + . ^r^^ll. i3™-i(A^), (49) 



where 

det £„(A^) = {dl -eA^c^x^)(6^ - e A^ x^)/x^. (50) 
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Substituting Am-i from this equation with m replaced by m — 1 into PS|) . we obtain a 
three-term recursion for Bm{^^)'- 

r. 



\ m—1 

N 



where we abbreviated 

'"m = '^m "~ ^m- (52) 

To define i3m(A^) by (jSH) we have to provide the initial values 
i3i(A^) = -eAV; i3o(A^) = 0. 

2.4- Fermat curve points appearing in the construction of the eigenvectors of Bn{X) 

As we have seen in the case of the two-site chain, the formulas (|38p . (j42p for the 
eigenvectors are given in terms of the points p2,o and p2 on the Fermat curve. The 
coordinates of these points are fixed by the values of amplitudes r2,o, ^"2,15 see In 
the n-site case four types of such points will appear: 

Pm {i-m-,ym) 1 Pm,s {_Xm,si Um.s) : Pm,s (5m,s ? ^m-.s) i Pm',s' m' ,s' V m' ,s'^ ' l^*^) 

The coordinates of these points are expressed in the terms of amplitudes rm,si m = 1, . . . , 
n, s = 0, . . . , m — 1 (defined as some solutions of equations (jTTj) . m = 1, . . . , n) by 

'^m' s' '^m,s / fm' ,s' ^ •^m,s (^m^m^m,s /brm ■^m.s dm/ (^^mCm^m,s) j S,S ^1. (^4) 

The corresponding y^fg,, ym,s are defined by the only condition on p^''^^,, pm,s to belong 
to the Fermat curve. The coordinates ym~i.h I < I < m — 2 are defined by 



m—2 m—l,l T-rm— 1 _ m,k 



m-l,l T-rm— i rri,« 

nl/m— l,s 1 lfc=l ym— 1,Z 17 1 n /rr\ 
T- 5 TT = 1, / = l,...,m — 2, (55) 
m—l,s T-rm— 3 m—l,l ' ' ' > V / 

s^l ym-l,l lls=l ym-2,s 

where 

The coordinates of the points Pmfi and pm are defined by 

(57) 

•^m^m '"m; yrrJ'm ^mdm^m—l/^m ■ (5^) 

Formulas are result from the solution of the eigenvector problem ((211), see 

Section 3. 
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The condition on the points Pm-i,i (1 < ^ < "^i — 2), p^.o and Pm defined by (jHS)), 
(jSZj), to belong to the Fermat curve gives 

m— 1 

^JV N TJ (^N _ „JV N _ N (,N _ „N N N n 
' m-l-^m-l' m,0 J_J_Vm-l,/ ' m,fcJ ~ 'm'm-2,0 I'^m-l "m-l ' m-l,«y' ^ 

^ ('"m-l,«Cm-l^m-l ~ '^m-l) JJ^ (^m-l,« ~" ''^m-2,s)5 / = 1, 2, . . . , m — 2, (59) 

s=l 

^iV _ /V ^A^^Af I j.iVj,JV . . . j,JV N 

fN ^ N N _ N ^ r'^<hNjN~N iN /^.x 

m m,0 m,l m,m— 1 m ' m m m—1/ m ' \^ J 

In order to relate these relations to the recurrent formulas of the classical BBS model 
1)511) we observe that the relations (jfiOj) (resp. ()fil)l ) follow from the relations obtained by 
the consideration of the highest (resp. lowest) terms in A in ()51|) starting from m = 2. 
Then, fixing in ()5T| successively at the m — 2 non- vanishing zeros of Bm-i, i-e. 
putting A^ = er^_ii, / = 1, . . . , m — 2, we obtain (jSHI)- Thus the points Pm-i,i, Pm,o 
and Pm defined by (jiS)), (jSTj), ()3H|) belong to the Fermat curve automatically. 

At the end of this section we would like to mention that the amplitudes rm,s can be 
found directly (i.e. not using the results from the previous subsection) from the relations 
fjKKjl . (jKTjl . (jK^ considered as equations with respect to rm,s and the coordinates of the 
Fermat curve points ()53p . These equations can be solved recursively starting from 
m = 2 and taking A^-th powers of these relations (see ()3I?|) . (jUUj) . (jHH))- Then the 
explicit formulas for the eigenvectors from the next section allow to obtain the Tarasov 
Proposition 1.5 in as a corollary. 



3. Inductive proof of the general solution of the auxihary problem 

Recall from p6|l that the vector \E'p^ (, := ipj^^a ^ Vi is eigenvector for i?i(A): 

5i(A)^,,„ = Ari,o^-''^'° ^Pi,o, 

and recall from (j^ . (jHUj) that the eigenvectors \l'p^ of -B„(A) were labeled by the vector 
Pn = (Pn,o, • • • ,Pn,n-i) £ C^nT- Let US further define: 

Pn = Efc=0 P",A:; Pn = (Pn,l, • • • , Pn,n-l) G (Zjv)"~^ (62) 

p^'^ denotes the vector p„ in which p„ ^ is replaced by pn,k ± 1, i-e. 

P^^ = (Pn,0, • • • , Pn,k ± 1, ■ ■ ■ , Pn,n-l), = 0, 1, . . . , n. 

The following Theorem ^ gives a procedure to obtain the eigenvectors \E'p„ G V''"\ 
n > 2, of -Bri(A) from eigenvectors ^E'p„_i € V''"^^'' of i?,„_i(A) and single site vectors 
ip^J^ G V„ defined by (jHS))- We start from \E'pi o- result of the first step of the 
induction we obtain the two-site result ()HHj) with P ^ .pH j) . 

The following theorem is valid provided ^ 0, the polynomials i3m(A^)/A^, 
m = 2, . . . ,n, have nonzero simple zeros and det 7^(er^^) 7^ (cf. the definition of the 
^-representation in P^). 
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Theorem 1 The vector 

= E Q(Pn-l>Pn|Pn)^p„_.®^i^ (63) 

P„_lG{Z^)"-l 

1 ,,(Pn-Pn-\){pn-Pn,0) 

^^^^^ Q{Pn^uPn\Pn) = J ^ X 

^ nr=l' mil ^p^!:Jpn-U - Pn,.) Wp^_,^X-Pn-u) ^g^^ 

n'^.'"=l (Pn-l.i - Pn-l,z) 7LT ^P^-i,XPn-l,l) 

is eigenvector of Bn{X): 

n.-l 

5„(A)^p„ = Ar„,o^-''-" n + ^n.-^^"'"'') ^Pn- (65) 

fc=i 

r/ie Fermat curve points p„, p„,i, p^j'^j and rn,k, entering (f^^P ^^^^^ related to the 
parameters of the model as,bs,Cs,ds, >i:s by equations |^j), / f33)) . |37| ), / f3^) . 
^ri(A) acts on \[^p^ as follows: 

n—l ^ X \ n—1 

An{x)^p^ = n ( 1 - ^ ) + ■ ■ ■ - ^"'^) + 



A — An.s \ A 



where ~ "-2 

¥'fc(p;) = F„(A„,,/c.) ni/^'i.s (67) 

s=l 

F„(A) = (6„ +cja„x„A) (Ac^ +dn/>in). (68) 
Corollary. In particular, at the n — l zeros Xn.k of the eigenvalue polynomial of Bn{X) 

K,k = -r„,fc^-^"■^ k = 1,. . . ,n-l (69) 

i/ie operator acts as shift operator for the k-th index of \l'p,j ; 

An {\n,k) *p„ = Mp'n) • (70) 

Further, the term in / I6'6|) of highest degree in X gives: V„ = V1V2 . . . v„ is a shift 
operator for the zeroth index of'^p^: 

V„ vl/^„ = vl/^+o . (71) 

Proof. We shall prove the Theorem Q by induction, showing that if it is valid for 
n — l site eigenvectors, then it follows for n site eigenvectors. Namely, we assume the 
correctness of the following formulas 

n-2 

Bn-i{X)^p^_, = Ar„_i,o^-^"-i'« J] (A - A„_i,0 (72) 

1=1 
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~ -^n-l,s \ 

^n-l,l ~ ^n-l,s I ^n-l,l 



1=1 \s^l 

n-2 , , X n-2 



+ n(l-T )^p„_, + Axi---x„_in(^-^"-M)-^p:°,' (73) 

s=l ^ n l,s/ ^^-^ 

where \n-\,i = —rn-i,i^~''"~^''' and the formulas for ipi{p'^_i) given by (jHTj) with n 
replaced by n — 1. 

Formula dHHI) for 5„(A)*p„: 

In order to prove the eigenvalue formula we use the following relation 

(A) f Aa„c„ + ^v„^ (74) 



which follows directly from (jHl) and pOj) . We apply its left-hand side to the left-hand side 
of (jU^ and its right-hand side to the right-hand side of (jHSl)- On the right, we use (f72j) . 
()73|1 ■ ()34|) ■ ()35|1 . According to ()73|1 . introduces shifts in the indices Pn-i 
while the second term involving v„ shift the index of ip'^^J . Since we are looking for an 
eigenstate, by shifting the summation indices we restore the original indices. However, 
this leaves a change in the matrix Q{pn-ii Pn\pn) ■ Now the difference equation 
for the tfp(7) functions appearing in Q{pn-ii Pn\Pn) is used, producing several factors 
under the summation which together we call R : 



5n(A)^p„= Q{pn-l,Pn\pn)R^p^.,®^^ 



in) 

Pn ' 

P„_ie(z^)"-i 

After some calculation we obtain 

R = ( TT A - \n-l,s I A ipil^p'"'- ) Pn\Pn) _^ 

s= 1 V An-l,.y ^^-^ (5(pn-l,Pn|pn) J 

n-2 



n., , .f. ^ Mn Q{pn-l,Pn - MPn) \ 
(A - A„_i,J Aa„c„ H — . — — 

V Q{pn-1, PnlPn) J 



+Ar„_i,ou; 
We have to show that 

n-l 

i? = A r,,ocu-^"-» JJ ( A - A„,fc) ; A,,^ = -r„,fca;-''".'= . (75) 
fc=i 

This will prove (jHSI)- Using the definitions of x^.s and Xm,s , we can rewrite 
for the argument A = Xn.k/^^ as follows: 

i^n(An,fc/o;) = \n,kbnCnUJ'^ (l " X„,fca;"''"''= ) (l - . (76) 
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Taking into account the expression for Q{Pn-i, Pn\pn), the definition for ^^(7) and 
the relations we obtain: 

0(n'^ n \n ) w^,k (p„-i,i - /0„,fc - 1) y, f^n + 

Qi.Pn-1, Pn\Pn) fjr ^ (p„_l,i - P„,fe) Wp^_^ X~ Pn-l,l) 



X 



"^Pn-l - 1) J7 V (Pn-l,s - Pn-l,Z + 1) i.»(Pn-l,/ " Pn-l,s " 1) 

' S^L \ ^n — l,s " n — l,l 



Q{pn~l, Pn\Pn) '^PnoiPnfi " Pn-l,o) ><1><2 " " " X„_ir„A;-/'" 

Q{pn-l,Pn - l|Pn) ^ ^,p„_i-p„ ^^Pn (Pn " Pn " 1) ^ _ f„^~^" - r^^'^" 

Q(Pn-l,Pn|Pn) ^^p„(Pn-Pn) ^nC^n f„_it^-^"-i 

Substituting these expressions into R gives 

(A - Xn^i,i) ■ Xa^Cn + ^ . 

After appropriate cancellations this becomes 

n-2 / \ \ \ \2 

« = E n /"V At'-"."--'"'" n(^..-.' - ^.a) + 

n— 2 n— 2 ^ , \ 

+ A2rr(A-A„_i,,)-r„.oa;-''-° + A f^o;"'^" TT ( 1 - ) . (77) 

To prove (f73j) we note that the coefficients at A" in both expressions (f73j) and (f77j) are 
rnflU!~''"'° and coefficients at A also coincide being rn^~^". Therefore the difference of 
these two expressions has the form A^P(A) where -P(A) is a polynomial of degree n — 3. 
Using the explicit expressions (fTHjl and (fTTjl we convince ourselves that P(A„_ij) = 
and therefore -P(A) = 0. This completes the proof that \E'p,j defined by (jUHj) . (jMj) is 
eigenvector of -B„(A) with eigenvalue (fT5|) . 

Formula ^ for A„(A„,fc)*p„: 

Next we show the validity of (|70|) . ()67|1 . We will need the relation 

u-\a^ - 6„v„,)A„(A) = (1 + Ax„cj-V„)5„(A)/A - v„ F„(A/cu) 5„_i(A)/A (78) 
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which can be obtained by ehminating between (I74j) and 

A„{\) = (1 + Ax„v„)A„_i(A) + u„(c„ - rf„v„)5„_i(A). (79) 
Let us apply (f7H|) to \E'p„ for A = — r„^fcCo'~''" *, i.e. at the zeros of Bn{\). This gives 

U.„;^(a„ - 6„V„)y4„(A„,fc)^p„ = - Fn{Xn,k/uj)/\n,k X 

X XI Q(Pn-l,Pn|Pn) 5n-l(An,fe)^p„_l (80) 

P„_ie(z^)"-i 

Pn GZ^ 

From (ff^ we know how to apply Bn-i to \I/p„_i: 

5n-i(An,fe)^p„_i = A„,fcr„_i,oa;-^"-i'0 JJ + r„_i,,a;-^"-i-) 



, (Pn-l,s - Pn,fe - 1)\ 

Using (jMjl we find the action of the inverse of the operator u^^(a„ — 6„v„) on ^/'p^'' : 

(u;i(a„ - 6„v„))-VJ? = r„-^c.^"<). (82) 
Shifting the summation over p„ in (j8(J|) and then applying (j82j) we obtain 

^n(A„,fc)^p„ = - ?^n^-^".(A„,fc/u;)/A„,fc X 

X Yl Q(Pn-l'Pn-l|Pn)^''"^n-l(A„,fc)^p„_, (83) 

p„_lG(Ziv)"-l 
PnGZjY 

Finally, using (|HT| and observing that 

Q(p„_i,p„-iip„)^^"-^"- n , — = ^-^^"+^-0 Q(p„_i,p„ip+'=) 

_i n,fc (p„_l,s- Pn,fc) 

we come to (jSJ). 
Formula (jTH) for V„^^p„: 

Using V„_i^p„_, = ^p+o^ and v^V'J,"^ = ^^^^^-^ we have 



V„_iV„^p,^ = X <5(Pn-l,Pn|Pn)^p+0^ ® 



(n) 

Pn + 1 



P„_ie(Ziv)"-i 

Pn GZjY 



X g(p;°i,Pn - l|Pn)^p„_, ® <\ (84) 



P„_ie(Ziv)"-i 

PnGZjY 



where in the second line we have shifted the summation variables pn-ifl and p„. Now 
considering the explicit form (jMj) for Q{pn-i, Pn\Pn) we read off that 

Q{p-'_„Pn - l\Pn) = Q{Pn-l,Pn\pi') (85) 

which gives (jTlj) . 

Formula 1^ for A„(A)^p„: 
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The operator A„(A) is a polynomial in A of nth order. From ^ and (|1(J|) we immediately 
read off its the highest and lowest coefficients: 

A,(A) = 1 + ... + A"xiX2...x„V„. (86) 

Using (jyij) we know how these terms act on ^ and if in addition we use the action of 
An at the n — \ particular values given in (fTOj) . we can reconstruct the action of the 
whole polynomial A„(A) on \l/p„ uniquely. Comparing this to we see that formula 

is the one which satisfies all these data. Therefore by uniqueness formula (jMj) is 
one which we are looking for. 

This completes the proof of Theorem ^ 



4. Action of on the eigenstates of 

In order to obtain the action of Dn{X) on we use the notion of the quantum 
determinant det^ T„(A) of the monodromy matrix. Since the rank of the matrix R{ujX, A) 
is 1, the definition of the quantum determinant is given by 

R{uX,X) T«(cuA)ri2)(A) = ri2)(A)r«(cuA)i?(cjA,A) =: detgT„(A) ■ R{ujX,\). (87) 
Explicitly we have 

detg T„(A) = Ar,{uj\)DrX\) - C„(a7A)5„(A). (88) 
Using pOj] and (jHTjl we obtain the factorization property of the quantum determinant 

detq Tn{X) = detg Li{X) ■ det, L2(A) ■ ■ ■ det, L„(A), 
For a single L-operator, using (j68|l . (j88|) gives detgLm(A) = VmFm(A) . So 

n 

An{u;\)Dn{\) - Cn{u:\)Bn{X) = V„ ■ J] F„(A). (89) 

m=l 

Acting by both sides of this identity on \E'p,^, fixing A = Xn,k (i-e. at the zeroes of the 
eigenvalue polynomial of -B„(A)) and using the inverse of (fTOj) with (jUTjl . we see that, 
very similar to An{\n,k), also Dn{\n,k) acts as a shift operator on '^p„, compare pTjl : 

Dn{Xn,k)^p„ = MPn) ^p+O'"^' ^'^(P'n) = ' ~^ ^nj^k^ H FMn,k). (90) 

^"-1 1L=1 yn-l,s m=l 

Note that Dn{\n,k) shifts pn,k in the opposite direction as A„(A„^fc) (see (lUTf) and (jZOI)) 
and Dn{Xn,k) also shifts p„^o- The shift in p„ q is due to the operator V„ at the right-hand 
side of dHHI)- Apart from the shifts just mentioned, applying the inverse of An{u!\) has 
cancelled in fj90|) the last factor m = n oi the quantum determinant ()89|). Analogously 
to (inni), using the particular values (plj) and reading off the coefficients of A° and A" 
directly from ()1()|1 . we obtain the following interpolation formula for D„(A)\l/p„: 

^n(A)v&p„ = g friT^^:^) ^Up'n)^pt-'^ + 
k=l \sj^k "'^ "'V '^''^ 

n— 1 ^ \ \ " L J " n— 1 

+ n 1 - ■ n ^ ■ + ^ n ■ - v.) ■ ^p.. (91) 

s=l ^ ^"'^^ m=l ^™ m=l s=l 
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5. Periodic model. Baxter equation and functional relations 

5.1. The Baxter equations 

After having determined the eigenvalues and eigenvectors of the auxiliary system, we 
now perform the first step of the program exposed in Subsection II. H[ i.e. the calculation 
of the eigenvalues and eigenvectors of the inhomogenous n-site periodic BBS chain model 
with the transfer matrix (fTT| . (fTl|l . Following the ideas of the papers [iHl El HH UHl IH] 
we are looking for eigenvectors of t„(A) as linear combinations of the eigenvectors "^p^ 
of the auxiliary system. 

It is convenient to go by Fourier transform in pnfi from \E'p^ to a basis of eigenvectors 
of V„ (and therefore of the Hamiltonians Hq and H„, see fll4|) . ()15p) 

^pX = Ep^.o^z^^-^-^-^^^P.; V„^,,^, = o;"^,,^,. (92) 

A shift of Pnfi in \E'p„ is replaced by a multiplication of ^E'p,p^ by powers of u. So from 
(!()() |1 and ^T\i . the action of t„(A) on ^E'p,pj^ becomes 



k=l \s^k 

+ {(i+-'n^)n(i 



, , X n-1 / . 



771=1 / s = l 

+ \(uJP f[ H^+Yl amCm I ~ ^"•^) \ ^ P'P'^^ ^^^^ 

\ m=l m=l / s=l J 

where we have taken into account that ipk{pn) Vk{p'n) independent of pn,o- Of 
course, since t„(A) commutes with V„, in there is no coupling between sectors of 
different p and we get separate equations for the different "charge" quantum numbers 
p which often will not be indicated explicitly. 

Let $p_E be eigenvector of t„(A) with eigenvalue 

t„(A| p, E) = + EiA + ■ • ■ + K_iA"-i + E„A" , (94) 

where E = {£'1, . . . , and from (fT3j) the values of £'0 and are 

n , , n n 

Eo =1+00^11^, E„ = n«™c^ n^-- (95) 

m=l ™ m=l m=l 

We are looking for $p_E to be of the form 

P'n 

From we get a difference equation for Q(p'„|p, E) with respect to variables 
which depends on A: 

t„(A|p,E)g(p;|p,E) = xi^^^-(/^n')Q(p';'lP,E)n^^^ + 

k=l ' Sytk ' ' 
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k=l s=ik 

n— 1 



771=1 

(?i n \ n— 1 ^ 

n ^™ + n ) n^^ - ^".-) \ Q^p'n\ 

m=l m=l / s=l J 

Substituting sequentially A = Xn,k, k = 1,2, ...,n — 1, we obtain a system of 
difference equations with respect to the variables p'„: 

tniKklP, E) Qip'J P, E) = f n ~ ] Mp'-J) QiP'n'l P, E) + 

+ fn^fr^^V'^Wr)WriP,E), A; = l,...,n-1. (98) 

In analogy to [iHl CHI HI we can decouple these difference equations using a Sklyanin's 
measure, namely, by introducing Q{pn \ P, E) defined as 

QM,,E)^^^!^ (99) 

lls,s'=i W n.s'{Pn,s — Pn,s') 



(si's') 



Rewriting (j^Hj) in terms of Q produces factors R± in both terms of the right hand side: 
UKAp, E) Q(p;|p,E)= I n ~ I a;-Vfc(Pn')i?-Q(Pn'|p,E) + 

+ f n ) ^'^V.(p'r) ^+ Q(P n'l P, E), A; = 1, . . . , n - 1.(100) 

where 

W n,fc(p„.s - Pn.k) WW-Mpn,k - Pn,s) 



J- W^n,k{pn,s - Pn,k - 1) Wp^^-s^{pn,k " Pn,s + 1) 



'^-l „n,fc 1 _ n,s p„^k-pn,.+l "~1 ,,n,fc x x x 

ni/n.s "''ri.A: ^ _ TT Wn.s '^n.s ^ ^n,s ~ ■^n,k 

1 n^fc _„ , 17^^^ ~ 11 o,"'* "\ ~\ \ 



and analogously -R_. We see that passing from Q to Q the brackets containing differences 
of terms Xn,i in (|100p are cancelled and so the equations decouple. This means that in 
terms of Q , the difference equations admit the separation of variables: 

n-l 

Qip'jE) = H UPn,k). (101) 

k=l 

Inserting the explicit expressions for (^^[p'^^) and Lpk{p'^^) we obtain Baxter type 
difference equations for the functions qk{pn,k)'- 

^n(An,fc|P, E) gfc(Pn,fc) = A+ ( A„_fe) ^fc (p„,fc + 1) + ^-{uj\n,k) qk{pn,k - ^) (102) 
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with , 

n— 1 

A+(A) = (.;7x.) (A/^)^-" n ^-(AH; A_(A) = Xfc (AM""^ F„,(AM; (103) 

m=l 

_ 71—1 n—2 

''"71,0 ''"n—1 /TT n.s / n,A;\ TT /in/i\ 
= ^ y (11 ^/n.fe/^/n.'s) 11 2/nll,.- (104) 

.=1 s=l 

In what follows we will mainly use t(A) instead t„(A|p, E). In fact the system of linear 
homogeneous equations ()102p with respect to qk{Pn,k), Pn,k ^ ^Af, is not completely 
defined. Since Ei, E2, . . . , En-i are unknown, the coefficients t(A„,fc) are also unknown. 
The requirement on the system of homogeneous equations ()102p for some fixed k, k = 1, 
2, . . . , n — 1, to have a nontrivial solution leads to the requirement that the matrix of 
coefficients must be degenerate. The latter gives a relation for the values Eq, Ei, . . . , 
En entering t(A). Taking all such relations corresponding to all = 1, 2, . . . , n — 1, and 
using the values of Eq and En given in (jHSI), at least in principle we can find the possible 
values of Ei, . . . , En~i- This fixes t(A). Then for every k, k = 1, 2, . . . , n — 1, we solve 
fll02j) to find qk{pn,k) for pn,k € Z^r (These difference equations have three terms and 
cannot be solved in terms of the functions Wp). This gives us finally Q{Pn\ P, E) and 
therefore the eigenvectors of the periodic BBS model: 

$p,E= ^-'-'"■'QiPn\P,^)'^P.- 

Pn = {pn,0, P'n) 

5.2. Role of the functional relations 

Now we will show that mentioned requirement on the systems of homogeneous equations 
fll02j) for all k to have a nontrivial solution is equivalent to functional relations EJ 
of the r(2)-model. We define r(o)(A) = 0, t^^\\) = 1, r(2)(A) = t{X) (see (jHl, dHSj)) and 

r(j'+i)(A) =t^^\uj'-^X)t^'\X)-uj'z{uj'~^X)t^^'^\X), j = 2, 3, . . . , iV, (105) 

where " 

^(A) =a;"''A+(A)A_(A) = J] F„(A/cu). (106) 

m=l 

The appearance of the monodromy determinant (jUHj) in the fusion relation is a direct 
consequence of the fusion procedure (see |22t l23]). 

The fusion hierarchy can be used to find eigenvalues of the transfer matrices in 
lattice integrable models. A key tools here is, in addition to ()105|) to demand a 
"truncation" identity which allows to express t^^\X) for some value j through r(*)(A) 
with i < j. A combination of the fusion hierarchy and truncation identity allows one 
to obtain an equation for t^'^\X) = t{X). This method was applied to many integrable 
models, in particular, to the RSOS models in and to the root of unity lattice 
vertex models in The functional relations for the r'-^-' -model for = 3 and the 

superintegrable case were first guessed in ^21 and have been solved to some extend in 
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The goal of the present section is to prove that the relations to determine the values 
El, . . . ,En-i entering t{X) also have the form of a truncation identity. We formulate 
this statement as follows: 

Theorem 2 The polynomial t^^'^^^X) satisfies the "truncation" identity 

^(iv+i)(A)_cu^2(A)r(^-i)(cuA) = A(A'^)+2^„(A^) (107) 

if and only if the system of homogeneous equations MU^) for all k has a nontrivial 
solution. 

Note, the classical polynomial A„(A^) + T'„(A^) corresponds to aq + cxq in (Tj. 

Proof. Let t{X) be a polynomial ()94|). ()95p such that the systems of homogeneous 
equations p02|l for all k have a nontrivial solution. We shall show that the polynomial 
P(A) = r(^+i)(A) - ujp z{\)t^^~^\uj\) at the left-hand side of is equal to 

-A„(A^) + ^^^(A^) With this aim we introduce the matrix 

r(A) = 1^ ^^'^(^^ 

Then it is easy to verify from ()105j) by induction that 

nu A) ■ ■ ■ r(u;A) r(A) - y -c.P.(A)rO-i)(..A) 
and we see that 

P(A) = tr r(cj^-^ A) ■ • ■ T{uj\) r(A) . (108) 

This relation shows the invariance of -P(A) under cyclic shifting A — > ujX. It means that 
in fact -P(A) depends only on A^. We denote V{\^) = -P(A). Thus we have to show 
that V{X^) = An{X^) +Vn{X^). The direct calculation gives that the coefficients of 
A° and A^" at both sides of this equation coincide. In order to calculate the coefficient 
in front of A° in the trace of the product of F-matrices (jlOSp one has to substitute 

where only the lowest terms in A in the matrix elements were kept. Therefore the lowest 
term in A in P{X) is 

-1 



tr 



Finding the eigenvalues of the matrix (|109p one can easily calculate the latter trace 
which is the lowest term 'P(O) and identify it with the lowest term 

n- l,N jN 

A(o)+p„(o)=i+n^ 

m=l 

of the polynomial ^^(A^) + Vn{X'^) calculated by means of relations and (jlUj) . The 
case of the coefficients in front of A^" can be treated analogously. 



Baxter-Bazhanov-Stroganov model: Separation of Variables and Baxter Equation 22 



Therefore to prove the Theorem |21 it remains to prove 

P{Kk) = A,(Ajfc) + 2^„(AJ,), = 1, 2, . . . , n - 1 , (110) 

where \n,k are given by (pTTjl and A^^^, = er^^i. are zeros of the polynomial (jUj). Let us 
fix some k and pn,k and denote the matrix of the coefficients of ()102j) with respect to 
the variables qkiPn,k), qkipn,k - 1), • • •, qkiPn,k - + 1) by M: 

( to -A^ ... -A+ \ 
-A+ ti -A^ ... 



M = -A 



V -^0 
where we abbreviated: = tiuo^X 



n,k ) 



■ ■ ■ ^N-l / 

Af = A±(uj^Xnk) ■ In order 



(111) 



to have a 



nontrivial solution, the matrix Ai must be degenerate. Let A4' be the matrix which 
has the same matrix elements as Ai except for Ai[ ^ = and A^^i = 0. Then, using 
the recursive definition ^QB) of t^^\\) and ^Q^, it is easy to show that the principal 
minor corresponding to the first j, j < N rows of the matrix Ai' gives T^-'^^\\n,k)- 
Calculating the determinant of the matrix pil|) and equating it to zero we obtain 

detM = T^^^'\\n,k) - C0fzi\n,ky~'\c0\n,k) 

- H A+(A„,fe^^)- n A-(A„,fc^^) = 0. (112) 



Further, 

n A.iXn,kCon = X^-irKt'^ detC^iXZk) 



Bn-l{Xn,k) 
' n '^n,k 



detCniKk): 



where we have used (ITnH|l . CT . dTUll . (IHill. (H7I) and 



Xk 



r5 (r^ 



-1 11 V n~l,s ' n,k 



n \ n,k/ g—\ n \ n,k/ 

Evaluating (jl^ at m = and A = \n,k so that Bn{Xnk) — ' finally we obtain 

n A_(A„,,c.^)=A(A^,). 

s e Z]v 

Substituting A = Xn,k into 

n 

det r„(A^) = A(A^) I)„(A^) - S„(A^) C„(A^) = det CM"") 

n rn=l 

= n n ^-(A^'"')= n ^(^^')= n (A+(Ac.'^)-A_(Aa;^)). (ns) 

m=l s e Z]v s e Zjv 

^^g^* n A4A.,...^)=I)„(A-,). 

Using (ITT^ we obtain (fTTTH) . 



s G Zjv 
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Conversely, if we have the polynomials t^^~^\X) and t^^^^\\) built from r(^)(A) = 
t[\) (see dnH), (in3)) by the recursion (jlOSj) and satisfying p07|) we get pi2|) at particular 
values of A. This means that the Baxter equations ()l()2j) have nontrivial solutions. 

This completes the proof of Theorem |21 



6. Periodic homogeneous BBS model for N = 2 

6.1. Solution of the Baxter equations 

In this section we consider in more detail the case of the N = 2 periodic homogeneous 
BBS model, where u = —1. By homogenous we mean that the parameters a, b, c, d 
and X are taken to be the same for all sites. As it was shown in ^H], for N = 2 and 
with arbitrary homogeneous parameters this model is a particular case ("free fermion 
point") of the generalized Ising model. 

We will find the eigenvalues tn{X\p, E) of the transfer-matrix t„(A) using a functional 
relation (see also [20], where a similar calculation is presented). We use the short 
notation t(A) for tn(A|p, E). From the previous section we have 

t(X) = 1 + (-1)^ + + ■ ■ ■ + ^„-iA"-^ + A"(a"c" + (-l)''x"). (114) 

Using for j = 2 and eliminating t^'^^ by p()7j) we get the functional relation 

t(A) t(-A) = (-l)^(^(A) + ^(-A)) + A(A') + V^X') (115) 

which we shall use to find t{X). Equivalently, we could have obtained pi5|) by 
multiplying together the two Baxter equations p02|) for A„_fc = ±r„ fc. 

Postponing for a moment the derivation (which will be supplied after ()126p ). let us 
anticipate that (|115|) can be rewritten as 

t{X)t{-X) = {-irl[{A{q)X' - C{q) +2iB{q)X), (116) 
q 

where 

A{q) = a c — 2xaccosg+x; B{q) = {ad — be) sin q; 
^/ X bd b'^d^ 

C g =1-2— cosg+^, 117 

X 

q is running over the set 7r(2s + 1 — p)/n , s = 0, 1, . . . , n — 1. Factorizing (|116p we get 

t(A)t(-A) = (-l)"J]A(g)(A -A,)(A +A_,) (118) 
q 

with 1 . „ 

\ = ]4(^( - ^m), D{q) = Aiq) C{q) - B{qf. (119) 

We fix the sign of ^jD{q) by the conditions 

^fD{q) = ^D{-q), g^0,7r; (120) 



^D{{]) = (x - ac)(l - M/x); ^ Din) = (x + ac)(l + M/x). (121) 
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In what follows we shall need the relations 

JJ A{q) = JJ(x - e'«ac)(x - e-'"^ac) = (a"c" + (-l)''x")2 , (122) 

J](/D(^ - iB{q)) = J](x - e"'ac){l - e'^bd/x). (123) 

q q 

The last relation can be obtained by grouping terms with opposite signs of q (modulo 
2-k) and using the definition of ^yD{q). Using p22|) we get 

t(A)t(-A) = (-l)"(a"c" + (-l)''x")2 - • + ^-'^) 

= (-l)"(a"c" + (-l)''x")2 J](A2 - A^), (124) 

where we made the change q —q in second product. From (jll4j) it follows that t{X) 
can be presented as 

n 

t(A) = (a"c'^ + (-l)''x'^)J](A-A,) 

s=l 

with zeroes A^. Therefore 

n 

t(A) t(-A) = (-l)"(a"c" + (-l)^x")2 JJ (A' - A^). 

Comparing with (jl24j) we obtain 

t{X) = {a-c- + (-l)''x") J](A ± A,), (125) 

9 

where the signs are not yet fixed. To fix these signs we compare the A-independent term 
in ()114|) with the corresponding term in p25p . The latter can be found using 

(a"c" + (-l)''x")nA, = (aV"+(-l)''x")-in(v^-ii?(g)) 

9 9 

= (a"c" + (-l)''x")-i JJ(x-e'''ac)(l-e'''6c//x) = (-1)^ + 6"c/"/x" , (126) 

9 

where we took into account p23|) . Therefore the number of minus signs in fll25|) must 
be even for the sector p = and odd for p = 1. Thus we have obtained 2" eigenvalues 
(2"'~^ each for both p = and p = 1). These eigenvalues provide the existence of 
nontrivial solutions of the system ()l()2j) of homogeneous equations. These solutions give 
the eigenvectors (|^. a basis in the space of states of the periodic BBS model for N = 2. 
We conclude this section supplying the derivation of ()116j) from ()115|1 : Using 

5+(A) := (6 + axA)(rf-cxA)/x, (5_(A) := 5+(-A) = (6 - axA)(rf + cxA)/x 

we easily verify the relations 61{X) = z{X), 5!! (A) = z{-X), 5+(A)5_(A) = (5(A2), 
where 2;(A) and 5(A^) are given by (jlOfij) and ()A.5|1 respectively. Taking into account 
An{X^) + Vn{X^) = tr {C{X^)f = xliX^) + x'liX^) and the relations 

r(A2) = tr/:(A') =x+(A=') + x_(A'), 5(A') = det ^(A^) = x+(A')x_(A'), 
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we can rewrite the functional relation (jlisp as 



t{\)t{-\) 



-l)P[z{X)+z{-X))+An{\^)+V^{\^) 

-l)P {51 + 51) +xl + x''_ 

-ly [xl + {-ir5l) {{x^/5^r + {-If) 

-lY\l{x^-e^^^){x^/5^-e^'^) 
-lY W ( e'^+ - t{\^) + e-'«5_ ) 



irn 



2 2 I 2\ \ 2 

a c + X )X — 



bd 2 
2 I A Hac] cosg 



—2iX{ad — be) sing}, 



(127) 



which confirms pifij) . 



6.2. Relation to the standard Ising model notations 

In the homogei 
(0) reduces to 



In the homogeneous N = 2 case we have cu = —1 and u^. ^ = u^, so the cyclic L-operator 



1 + Axvfc Aufe (a - 6vfc) 
Ujt(c— c/vfc) Xac + Vkbd/x 

Let us make the special choice of the parameters d = be/ a and A = 6/ (ax). Then 

1 Ukb/ X 



128) 



Lfc(A) = (l + Vfe6/a) , 

CUfc oc/x 
and the transfer-matrix is 



;i + Vfc6/a) I ^ j C 1 Ufe6/x 



(129) 



t„(A) = tr Li(A) L2(A) ■ ■ ■ L„(A) = J](l + ■ b/a) ■ J](l + u,_iu, • 6c/x). (130) 

k=l k=l 

Recall that due to the periodic boundary conditions u^+i = u^. Using 
exp{KiUk_iUk) = coshi^i(l+ u/j._iUfctanhi^i); exp(i^"2Vfc) = coshK2(l+ '^fctanh/r2)5 

and writing = cr| and = erf, it is easy to identify t„(A) with the standard Ising 
transfer-matrix: 

K2 \ exp I ^ Ki al_^ al \ ; tanhK^ = -; tanhi^i = — . 



7. Conclusion 



This paper is devoted to the solution of the eigenvalue and eigenvector problems for the 
finite-size inhomogenous periodic Baxter-Bazhanov-Stroganov quantum chain model. 
We use an approach which had been developed in full detail for the quantum Toda 
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chain in j^l CD] and in PQ] for the relativistic deformation of the Toda chain. This 
approach consists of two main steps: In order to find eigenvectors for the transfer matrix 
^n(A) + Dn{\) first we find the eigenvectors of the off-diagonal operator -B„(A) adapting 
the well known recurrent procedure described in PU] to our root-of-unity case. Then, 
using these eigenvectors, we construct the eigenvectors for the BBS transfer matrix and 
show that the coefficients of the decompositions of one set of eigenvectors in terms of 
the other set factorizes into a product of single variable functions, each satisfying the 
Baxter type equation. We show that the condition for these equations to have nontrivial 
solution is equivalent to the functional relations for the transfer matrix eigenvalues in 
the BBS or r*^^^ model. In case of = 2 the Baxter equation can be solved and as result 
we obtain the eigenstates of the transfer matrix of the generalized Ising model at the 
free fermion point. We shortly give the relation of the N = 2 BBS model parameters to 
the standard Ising model parametrization. 
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Appendix: Amplitudes rm,k in the homogeneous case 

The determination of the amplitudes r^^k for the inhomogenous BBS model had been 
reduced to solving Equation pTjl with (fFT|) of Section 12.31 Here we show that in the 
homogenous case this task simplifies to solving just one quadratic equation, using a 
trigonometric parametrization. 

In the homogeneous case we have 

am = a, bm = b, Cm = c, dm = d, Km = X, = r, CmiX^) = C{\^) (A.l) 

and 

C„(A^) P.(A^) ) - ^^^^ ■ ^^-2^ 
Using the fact that a 2 x 2 matrix M with eigenvalues /i+ and /i_ satisfies 

jYjm ^ /^+ jyj _ ^^+^^- ^^~^^+ ^ 

from the matrix C{\^) we obtain 

i3^(A^) = -eAV^i-^, 

— X— 

where x+{\^) and a;_(A^) are the eigenvalues of £(A^). These eigenvalues are the 
roots of the characteristic polynomial — t(A'^) x + 5(A^) = : 

x± = l(r± Vr2 -45), (A.3) 
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where, see (jSOI), 



r(A^) = tr /:(A^) = 1 + ^ - 6A^(x^ + a^c^), (A.4) 

(5(A^) = det /:(A^) = - eA^ a^) (rf^ - eA^ x^) . (A.5) 

Introducing the variable by = e"^ we find that roots of Bm correspond to 

roots (t)m,s of e'"^'^ = 1 (without = 0): 

4'm,s = '^'^s/m, s = 1, 2, . . . , m — 1. (A. 6) 

Now we need to find the exphcit relation between A^ and 0. We have 



r + Vr2-4(5 = e'"^ (r - Vr^ - 45) or r2 = 45cos^^. (A.7) 

Taking into account ()A.4|) and ()A.5|) . we consider ()A.7|) as a quadratic equation for A^: 
A2^(a2^c2^ + x^^ - 2a^c^x^ cos 0) + (fo^^rf^^ + x^^ - 26^rf^x^ cos 0) /x^^ 

- 2eA^ (^(a^ - &^)(c^ - ^^) + ^"^^^f "^"^ + - («^^^ + cos0^ = . (A.8) 

The solution A^(0) of this equation describes the relation between the variables A^ and 
0. Therefore we can translate the zeroes ()A.6|) of Bm{^^ {4>)) in terms of variable 0, to 
zeroes A^(0m_s). From (jlTjl we get 

rZ,s = eA^(0^,.), s = 1, 2, . . . , m - 1. (A.9) 

The value of q can be found recursively from (jHIH) using (jA.ljl : 

<o = a'^c^ + , r^, = r^ = a^-h^. (A.IO) 
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